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Abstract. Recent hadronic charm decay results from fixed-target experiments are presented. New 
measurements of the D'^ K^K^K^K^ branching ratio are shown as are recent resuhs from Dalitz 
plot fits to D+ K'K+%+, 7t+7t"7i:+ K'%+%\ K+%^%+ and D+ K+%'%+. These 

fits include measurements of the masses and widths of several light resonances as well as strong 
evidence for the existence of two light scalar particles, the TtTt resonance O and the Kn resonance K. 



INTRODUCTION 

Hadronic decays of charm are rich in information about QCD. For instance, the sup- 
pression of TT+TT^ relative to K proved the importance of final state 
interactions in charm decays. Also, hadronic decays give rise to the 2.5 x difference 
between the D+ and lifetimes. Hadronic decays can provide information on relative 
strengths of decay diagrams (spectator, exchange, annihilation, etc.). Spectator diagrams 
are believed to be responsible for the vast majority of the charm decay rate. In a spectator 
diagram, the charm quark decays while the other quark in the meson is a spectator. By 
contrast, exchange and annihilation diagrams require a connection between the charm 
quark and the other quark in the meson and are therefore suppressed. Determining the 
contributions of these diagrams is an interesting open question in charm physics. 

More recently, charm has been used to investigate the light resonances which are 
products of charm decays. Although very high statistics scattering experiments have 
been performed for many years to investigate these resonances, many parameters are 
still virtually unknown. Charm offers a unique way to investigate these resonances by 
nature of its low background and well defined initial state (pseudoscalar meson). 

Experimentally, hadronic decays can be investigated by comparing branching ratios 
and by analyzing the resonant structure of multibody decays. The results presented here 
come from the Fermilab experiments E791 and FOCUS. E791 (FOCUS) ran in 1991 
(1996-7) with a 500 GeV/c^ TT" beam (180 GeV photon beam) on five (four) targets. 



BRANCHING RATIO MEASUREMENTS 

Calculations of charm quark decay rates via the weak interaction have been possible for 
many years. Unfortunately, only charm hadron decay rates (which are affected by the 
strong force) can be measured by experiments. The strong force is even more intimately 
involved when the charm particle decays into hadrons. Thus, deviations from the naive 
weak prediction for a given decay can provide insight into the nature of the strong force. 
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FIGURE 1. E791 and FOCUS signals from D" K-K+K-%+ and D" K-%+%-%+ decays. 

K~K^K~n^ branching ratio 

The decays K+K and K'Ti+Ti^n+ are both Cabibbo favored de- 

cays. A Cabibbo-favored hadronic D'^{cu) decay produces one s quark from the charm, 
one u spectator quark, and a u and d quark from the virtual W. If both decay modes 
were entirely non-resonant, forming the K^Tt^n^Tt^ final state would require pop- 
ping uiidd from the vacuum while K K^K would require uuss and the branch- 
ing ratio between — > K K^K and — > K^n'^n^n'^ could be used to deter- 
mine the ss suppression relative to dd. Multi-body charm decays generally proceed 
through resonances, however, which complicates the issue. Even in this case, one can 
note that the ^^71+71^71+ decay can occur through resonances using only the four 
quarks from the decay, e.g. K*{S92)^{sd) p{nO)^{uu), while the K+K %+ 

decay requires either an ss pair from the vacuum or a final state interaction which cou- 
ples rni to KK. The recent E791 result contains some of these speculations Sig- 
nals for these two modes from E791 and FOCUS are shown in Fig. |l]. E791 finds 

^^{ jD^K-%^-l+ ) = 0-0054 ±0.0016 ±0.0008, significantly higher than the E687 
measurement of 0.0028 ± 0.0007 ± 0.0001 0]. The preliminary FOCUS resuk is 
much closer to the E687 result at 0.00306 ± 0.00047 (statistical error only). 

DALITZ PLOT ANALYSES 

Multibody decays of charm particles can occur through various strong resonances which 
can interfere with each other. In a three-body decay, a "Dalitz" plot can be constructed 
to show the effect of these resonances and their interferences by plotting the squared 
invariant mass of two combinations of the final state particles against each other. In the 
absence of interference, how a resonance appears on the Dalitz will depend on its mass 
and width (a relativistic Breit-Wigner) as well as on its spin (Legendre polynomials). 
Interference effects can significantly alter these shapes. Fitting a Dalitz plot to a fully 
coherent sum of resonances allows one to extract information about how much each 
resonance contributes and how each resonance interferes with other resonances. By 
performing a coherent Dalitz plot analysis one can extract information about weak 
decays and the effects of the strong force on the weak decays. 
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FIGURE 2. Preliminary K /r+7t+ FOCUS Dalitz plot and projections. Projections show the data 
background, data signal+background, and the fitted result. 



TABLE 1. Preliminary FOCUS Dalitz plot fit results 
for the decay K^K^k^. (Statistical errors only.) 



Decay mode 


Fraction (%) 


Phase (°) 


K*{892)K+ 


22.0±1.1 


(fixed) 


ao{980)n+ 


27.8±4.8 


146 ±5 


(j)(1020)7i:+ 


27.8 ±0.9 


244 ±6 


/2(1270)7t+ 


0.7 ±0.2 


12±7 


/o( 1370)71+ 


5.9±1.2 


60±6 


K*{U10)K+ 


8.8±1.9 


135±6 


K^{U30)K+ 


69.3±6.3 


63±4 


(^(1680)71+ 


1.5±0.5 


-70 ±9 


sum 


163.8 





D+ K-K+71+ decays 

FOCUS has a high statistics sample of the singly Cabibbo suppressed decay 
K K^n^. The Dalitz plot as well as projections along each axis for this mode are shown 
in Fig. ^ In the Dalitz plot and on the nP-{KK) projection, a very clear (|)(1020) can be 
seen. The ^*(892) across the Dalitz plot and along the irP'{Kn) projection is also clear. 
Both contributions show a cosG modulation of the amplitude due to the spin-1 nature of 
the resonances. Distortions due to interferences are also visible. Preliminary fit results 
from a fully coherent analysis of the data are tabulated in Table |l]. These results support 
the presence of significant amounts of ^* (892) and (|)( 1020) . The large contribution from 
the ^q(1430) partly explains the broad enhancement at high nP-{Kn). The existance of 
many resonances with very different phases indicates significant interferences, as does 
the fact that the fit fraction sum is much greater than 100%. These interferences explain 
the obvious distortions seen in the Dalitz plot. Work is currently underway to investigate 
direct CP violation by comparing the Dalitz plots of the Z)+ and D decays. While for 
two-body decays a simple branching ratio is sufficient, for multi-body modes a Dalitz 
plot analysis is necessary to extract all of the information on direct CP violation. 



FIGURE 3. Diagram of annihilation (left) and resonance (right) contributions to D+ — > 7t+Jt Jt+. 




FIGURE 4. Preliminary Df 7t+7t 7t+ FOCUS Dalitz plot and projections. Projections show the data 
background, data signal+background, and the fitted result. 



Dj 71+71 71+ decays 

Although D+^TT+TC 71+ is Cabibbo favored, it can only occur via a spectator diagram 
if it uses a resonance which couples to both KK and nn or via final state interactions. It 
can also proceed via an annihilation diagram. These possibilities are sketched in Fig. |^. 
Since p(770) does not couple to KK, any p(770) would indicate an annihilation diagram 
contribution or final state interactions. A resonance known to couple to KK and nn is the 
/o(980). This mysterious state has been proposed as a four-quark state and KK molecule 
among other things. The presence of the ao(980) further complicates the understanding 
of this unique resonance. The /o(980) mass is below the KK threshold but is broad 
enough to have a significant branching fraction to KK even with limited phase space. 

The FOCUS Dalitz plot and projections for 7i;+7i;^7C+ are shown in Fig. ^. 

Clear /o(980) bands are visible in the Dalitz plot and projections. A concentration at 
m^(7r+7r^) ~ 2 GeV/c^ is also visible. The results of a preliminary fit to this distribution 
as well as E791 results [|3|] are shown in Table 0. Both results clearly show /o(980) 
dominance and no significant p(770). Thus, this decay proceeds almost entirely through 
resonance modes with no evidence of an annihilation diagram contribution. 

Using data, E791 and FOCUS have made measurements of the poorly measured 
scalars contributing to the D+ 7r+7r^7r+ decay. Although the PDG fists /o(980) 
mass and width measurements from 1973 to the present, with 13 measurements in 1999 
alone, there is still no consensus on either value. The 1=0, J=0 states betweeen 1200 and 
1500 MeV/c^ are even more murky. 

E791 finds amass and width of the /o( 1370) of 1434± 18±9 MeV/c^ and 172±32± 
6 MeV/c^, respectively. FOCUS uses 5o( 1475) for a scalar around 1475 MeV/c , as seen 
in an E687 analysis. The preliminary mass and width are found to be 1473 ± 8 MeV/c^ 
and 1 12 ± 17 MeV/c^, respectively; quite comparable to the E791 result for /o( 1370) . 



TABLE 2. E791 and preliminary FOCUS Dalitz plot fit results for the decay 
Tt+TfTt"'". Errors on FOCUS results are statistical only. 



E791 FOCUS (preliminary) 



Decay mode 


Fraction (%) 


Phase (°) 


Fraction (%) 


Phase (°) 


/o (980)71+ 


56.5±4.3±4.7 


(fixed) 


94.4 ±2.5 


(fixed) 


non-resonant 


0.5±1.4±1.7 


181±94±51 


25.5 ±4.4 


246 ±4 


p (770)71+ 


5.8±2.3±3.7 


109±24±5 






/2(1270)7t+ 


19.7±3.3±0.6 


133±13±28 


9.8±1.2 


140 ±6 


p(1450)7i+ 


4.4±2.1±0.2 


162 ±26 ±17 


4.1 ±0.7 


188 ±14 


5o(1475)7i:+ 






17.4 ±2.2 


250 ±4 


/o( 1370)71+ 


32.4±7.7±1.9 


198±19±27 






sum 


119.3 




151.2 





E791 finds slightly better fits using the WA92 coupled channel Breit-Wigner for- 
malism to describe the /o(980). A standard relativistic Breit-Wigner is proportional 
to l/(m^ — + imrTr{m^)) for a resonance r with a mass and width of m,- and F,-, 
respectively and at a two-body mass- squard of m^. In the /o(980) WA92 coupled 
channel formula, is replaced with -|- where T^ijrP-) = ^ti A/m^^^T^'-'wi 

and Ff(m^) = Sk{\J fn%n/^ ^ '^a:+ + \J '^^jr/ ^ ^ "^io j ■ This framework yields 

My„(980) = 977 ± 3 ± 2 MeV/c2, g^, = 0.09 ± 0.01 ± 0.01, and gK = 0.02 ± 0.04 ± 0.03 
from the E791 data. Fitting with a standard Breit-Wigner results in minor changes: 
MyQ(980) = 975 ± 3 ± 2 MeV/c2 and Fy,,(98o) = 44 ± 2 ± 2 MeV/c^. FOCUS finds 
the K-matrix framework works quite well in dealing with the coupled channel 
nature of the /o(980) In this framework, transformed variables are used: 

ml = + (yKK/yTZTt) {\9KK{mr)\/pTtK{mr))mrrr and Fq = m^F7(mopjt7i(m^)7^j^) 
where and pKK are phase space terms and Jnn and Jkk are coupling constants 
normalized to Hn'^yKK = 1- In this framework, the preliminary FOCUS fits return 
M/o(980) = 963 ± 6 MeV/c2, F^„(98o) = 297 ± 92 MeV/c2, and f^Jy^^ = 2.09 ± 0.53 
which translates to My„(98o) = 982 ± 30 MeV/c^ and Fyy(98o) = 89 ± 32 MeV/c^ for a 
standard Breit Wigner (errors are statistical only). 



D+ 71+71 71+ 

E791 has published results for a coherent Dalitz plot analysis of 7i;+;c 7C+ 
In their initial fit to the Dalitz plot using all known resonances, the fit quality was very 
poor with a confidence level of 10^^. By including a low mass scalar particle (a) the 
fit was significantly improved and yielded a confidence level of 75%, providing strong 
evidence for the elusive light scalar. The results from both fits are tabulated in Table |3|. 
From the fit, the o parameters were determined to be = 478_23 i 17 MeV/c and 
Fa = 324til ±21 MeV/c2. Many checks were made to validate the existence of the o in 
this decay. These checks included fitting with a vector and tensor state, and fitting with 
no phase variation. The fit with a phase-varying scalar particle was clearly preferred. 



TABLE 3. E791 Dalitz plot fit results for the decay D+ 71+71^71+ Param- 
eter errors from "without o" fits are statistical only. 





Without a- CI 


- 001 % 


With a: CL = 75% 


Decay mode 


Fraction (%) 


Phase (°) 


Fraction (%) 


Phase n 


p (770)71+ 


20.8 ±2.4 


(fixed) 


33 6±32±22 


ffixed^ 


non-resonant 


38.6±9.7 


150±12 


7.8 ±6.0 ±2.7 


57±20±6 


jQ[yoV)Tl 


/ .4 zt i .4 


1 ^0 1 1 


6.2 ±1.3 ±0.4 


165 ±11 ±3 


/2(1270)7t+ 


6.3 ±1.9 


103 ±16 


19.4±2.5±0.4 


57±8±3 


/o(1370)7t+ 


10.7±3.1 


143 ±10 


2.3 ±1.5 ±0.8 


105±18±1 


p(1450)7i+ 


22.6 ±3.7 


46 ±15 


0.7±0.7±0.3 


319±39±11 


07t+ 






46.3±9.0±2.1 


206±8±5 


sum 


106.4 




116.3 










^71+ 





The Cabibbo favored decay K^n'^n'^ provides a very high statistics mode in 

which to study charm decays. Previous analyses of this decay ^ have identified two 
mysteries in this decay. The first mystery is why there is a dominant non-resonant con- 
tribution; unique in charm decays. The second mystery is why a good fit to this Dalitz 
plot seems to be impossible to achieve. The E791 data sample of 15,090 events (94% 
signal) can be used to shed light on these mysteries. Fitting the data using all known res- 
onances results in a large non-resonant contribution and a very poor fit (confidence level 
of 10^^^), similar to past attempts. Given the evidence for the a in the tc+tc^tc^ 
mode, an additional scalar (k) was added to the — > Z^7r+7r+ fit. This provides a 
dramatic reduction in the non-resonant contribution (from 91% to 13%) and a much im- 
proved fit (confidence level of 95%). The preliminary results of both fits are shown in 
Table ^. The preliminary K parameters are found to be Mk = 797 ± 19 ± 42 MeV/c^ and 
Fk = 410 ± 43 ± 85 MeV/c^. Attempts to explain the data in other ways (e.g. includ- 
ing a vector state, tensor state, non-phase- varying state, structured non-resonant contri- 
bution) have been inadequate. Preliminary meausrements of the ^^(1430) parameters 
(^i^o*(i430) = 1459 ± 7 ± 6 MeV/c2 and Fj^*(i43o) = 175 ± 12 d= 12 MeV/c^) have also 
been extracted in this analysis [1T0|]. 

TABLE4. Preliminary E791 Dalitz plot fit results for the decay D+ ^/:"7t+7i:+. 
Parameter errors from "without k" fits are statistical only. 

Without k: CL = 10^ ' With k: CL = 95 % 



Decay mode 


Fraction (%) 


Phase (°) 


Fraction (%) 


Phase (°) 


NR 


90.9±2.6 


(fixed) 


13.0±5.8±2.6 


349 ± 14 ± 8 


F (892)71+ 
r;( 1430)71+ 
r2( 1430)71+ 
F (1680)71+ 


13.8±0.5 


54±2 


12.3±1.0±0.9 


(fixed) 


30.6±1.6 


54±2 


12.5 ±1.4 ±0.4 


48±7±10 


0.4±0.1 


33±8 


0.5±0.1±0.2 


306±8±6 


3.2±0.3 


66±3 


2.5 ±0.7 ±0.2 


28 ±13 ±15 


K7t+ 






47.8± 12.1 ±3.7 


187±8±17 


sum 


138.9 




88.6 






FIGURE 5. FOCUS Dalitz plots and projections for D+ K+%-%+ and 
show the data background, data signal+background, and the fitted result. 
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TABLE 5. Preliminary FOCUS Dahtz plot fit results for 
D^—^K^K~K~^. Statistical errors only. 



■ K^K 71+ and 



Decay mode 



D+ K+%-%+ 
Fraction (%) Phase (°) 



Fraction (%) Pliase (°) 



p (770)7*:+ 


51 ±10 


(fixed) 


40 ±4 


(fixed) 


non-resonant 


9±5 


-6±16 


18±4 


34±7 


(892)71+ 


43±7 


208 ± 16 


22 ±3 


163 ±7 


/o(980)/:+ 


9±5 


73±31 






/2(1270)/:+ 






2±1 


33±21 


K*{U\Q)%+ 


12±8 


133±23 


14±5 


-10±7 


/sTq* (1430)71+ 






14±6 


68±7 


/s:|( 1430)71+ 


6±3 


48 ±27 






p(1450)/:+ 


10±5 


247 ± 15 


8±2 


219±14 


K*{\6%0)%+ 


22 ±10 


2±20 






sum 


162 




118 





FOCUS has obtained preliminary results from Dalitz plot fits to the doubly Cabibbo 
suppressed decay D+ — ^ K^n^Tt^ and the singly Cabibbo suppressed decay Df 
K^n^Tt^ . The Dalitz plots and projections are shown in Fig. |5| and the fit results in 
Table |5|. This is the first fit to the K'^n^Tt^ Dalitz plot. Both fit results indicate a 

rich resonance structure, dominated by p(770). 



TABLE 6. Fitted masses and widths for various scalars involved in charm decays. E791 K 
and ^"0(1430) results are preliminary. The FOCUS results are preliminary and errors shown 
are only statistical. 



E791 



FOCUS 



Resonance 



Mass (MeV/c2) F (MeV/c^) 



Mass (MeV/c2) F (MeV/c^) 



a 478+15 ±17 

K 797 ±19 ±42 

/o(980) 975 ±3 ±2 

/o(1370)/5o(1475) 1434±18±9 

K*{U3Q) 1459 ±7 ±6 



324+^^ ±21 

410±43±85 

44±2±2 

172±32±6 

175±12±12 



982 ±30 
1473 ±8 



89 ±32 
112±17 



CONCLUSION 



Hadronic decays of charm provides an environment to study many aspects of high energy 
physics including measuring the contributions of various Feynman diagrams, studying 
the effect of final state interactions, searching for CP violation, and measuring the mass 
and width of light resonances. The relatively large decay rate of — > n'^n^n'^ decay 
has been found to be due to resonances which couple simultaneously to KK and nn 
rather than annihilation diagrams. Strong evidence for, and precise measurements of, 
two particles which have existed on the fringe for many years (o and k) are presented. 
Using the clean charm environment to measure light resonance parameters is a new and 
interesting use of charm hadronic decays. Table ^ summarizes all of the measured light 
resonance values described in these proceedings. The future holds the prospect of even 
more precise measurements of light resonances, a search for direct CP violation from 
Z)+ — > K K^n^ decays, and much more information from many decay modes. 
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